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Cn , Abstract. A factorization of a permutation into transpositions is called "primitive" if its factors are weakly ordered. 

We discuss the problem of enumerating primitive factorizations of permutations, and its place in the hierarchy of pre- 
viously studied factorization problems. Several formulas enumerating minimal primitive and possibly non-minimal 
primitive factorizations are presented, and interesting connections with Jucys-Murphy elements, symmetric group 
characters, and matrix models are described. 

j£^ • Resume. Une factorisation en transpositions d'une permutation est dite "primitive" si ses facteurs sont ordonnes. 

j^ ' Nous discutons du probleme de I'enumeration des factorisations primitives de permutations, et de sa place dans 

la hierarchie des problemes de factorisation precedemment etudies. Nous presentons plusieurs formules enumerant 
certaines classes de factorisations primitives, et nous soulignons des connexions interessantes avec les elements Jucys- 
Murphy, les caracteres des groupes symetriques, et les modeles de matrices. 

►^ , Keywords: Primitive factorizations, Jucys-Murphy elements, matrix integrals. 
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1 Introduction 

lO ' The problem of counting the number of ways in which a given permutation can be factored into a given 

^P , number of transpositions is of perennial interest in algebraic combinatorics. Usually, one considers this 

problem in the presence of constraints on the factors, e.g. that they should be transpositions of a cer- 
tain type, or should collectively generate a certain group, etc. By varying these constraints, one obtains 
enumeration problems which enjoy surprising connections with other branches of mathematics. 
K> , The earliest enumerative study of transposition factorizations was carried out by Hurwitz [|l5|] in the 

^ ' nineteenth century. Motivated by a problem from enumerative algebraic geometry, namely the counting 

of almost simple ramified covers of the sphere by other Riemann surfaces, Hurwitz published an explicit 
formula for the number of minimal transitive factorizations of an arbitrary permutation into transpositions. 
There are two constraints in the Hurwitz factorization problem: "minimality," which requires that the 
number of transpositions used should be as small as possible, and "transitivity," which requires that the 
factors should act transitively on the points {!,...,«}. Hurwitz's formula for the number of minimal 
transitive factorizations of a permutation vr G S{n) of cycle type ji — {jjli , . . . , ^e) h n is 
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As a particularly beautiful special case, Hurwitz's formula yields that the number of factorizations of a 
full cycle in S{n) into n — 1 transpositions is 71"^^, the number of trees on n labelled vertices. This case 
of Hurwitz's formula was independently rediscovered and popularized by Denes [pi] . The general Hurwitz 
formula was independently rediscovered by Goulden and Jackson JlCJ], to whom the first rigorous proof 
is due. 

A key feature of the Hurwitz factorization problem is centrality: the number of minimal transitive 
factorizations of tt depends only on the cycle type of tt. This remains true for transitive factorizations 
of arbitrary length. A different choice of constraints leading to a non-central factorization problem was 
considered by Stanley p^, who initiated the study of what he termed reduced decompositions. These are 
minimal factorizations in which the transpositions allowed to be used as factors are the Coxeter generators 
(s, s + 1). Reduced decompositions are also referred to as sorting networks because of their relation to 
the bubblesort algorithm familiar to computer scientists. They could also be called minimal Coxeter 
factorizations. The enumeration of reduced decompositions is complicated by its non-centrality, and has 
spawned its own extensive literature, see |S] for a beautiful introduction. The asymptotic behaviour of 
random reduced decompositions is the subject of an intriguing set of conjectures due to Angel et al. [|l|]. 

The Coxeter factorization problem naturally fits into a wider class of constrained factorization prob- 
lems, in which the factors are chosen from a specified set of transpositions which generate S{n). A 
second example from from this class was considered by Pak [E7|, who initiated the study of star factor- 
izations. These are factorizations in which the transpositions allowed to be used as factors have the form 
(1*). For example, the unique minimal star factorization of (123) G S'(3) is (123) — (13)(12). Thanks 
to recent work of Irving and Rattan [p^, Goulden and Jackson [pi]], and Feray [^, the combinatorics of 
star factorizations is now completely understood. 

Recently, Gewurz and Merola [0] posed the problem of enumerating transposition factorizations under 
a constraint on the order of the factors. A factorization 

■n = {.9i,ti) ...[sk.tk) [si<ti) (2) 

of TT e S{n) into a product of k transpositions is called primitive if 

2 < ti < ■ • • < t/c < n, (3) 

i.e. if its factors appear in weakly increasing order with respect to the larger element in each. For example, 
(123) e S{3) can be factored into a product of two transpositions in three ways, 

(123) = (12)(23) = (23)(13) = (13)(12), (4) 

but only the first two of these factorizations are primitive. Gewurz and Merola obtain the interesting result 
that the number of primitive factorizations of the cycle (12 . . .n) E S{n) into n — 1 transpositions (i.e. 
the number of minimal primitive factorizations) is the Catalan number 

1 /2n - 2\ 
Cat„„i = - . (5) 

n \n — I / 

This result should be considered in tandem with Hurwitz's n"^^-count of unrestricted factorizations of 

(12... n). 
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In this extended abstract prepared for FPSAC 2010, we will give an overview of the authors' ongo- 
ing work on the enumeration of primitive factorizations. First we study minimal primitive factorizations, 
and obtain an analogue of Hurwitz's formula (lib, i.e. an explicit formula which counts minimal prim- 
itive factorizations of an arbitrary permutation (Theorem n^ and Corollary n^ below). Then we present a 
link between the primitive factorization problem and Jucys-Murphy elements. This connection explains 
the centrality of the primitive factorization problem, and allows us to use character theory to enumerate 
primitive factorizations of a full cycle into any number of transpositions (Theorem below). Finally, we 
discuss a surprising connection between the primitive factorization problem and the theory of matrix mod- 
els: generating functions enumerating primitive factorizations may be expressed as integrals over groups 
of unitary matrices against the Haar measure (Theorem]^ below). It turns out that these integrals are of 
independent interest and have a long history in mathematical physics. 

2 Minimal Primitive Factorizations 

Any primitive factorization of tt into k transpositions has the form 

■n = (*2) . . . (*2) (*3) . . . (*3) . . . (*n) . . . (*n), (6) 



where (02, . . . , a„) is a weak {n — l)-part composition of k. We will say that the above factorization is 
of type A h fc if the frequencies a2, 03, . . . , a„ coincide with the parts of A after reordering. For example, 
there are three primitive factorizations of (1234) £ S{A) of type (2, 1), namely 

(1234) = (23)(13) (34) = (12) (34)(24) = (23) (34)(14) . (7) 

2112 12 

Let us now enumerate minimal primitive factorizations by type. Let (£(A:) denote the set of all weakly 
increasing sequences ii < ■ ■ ■ < ik of k positive integers such that ip > p for 1 < p < k — 1 and ik — k. 
It is not difficult to show that 

|(£(fc)|-Catfe. (8) 

Given a partition A h fc, one may introduce a refinement RC(A) of the Catalan number by declaring 
RC(A) to be the number of sequences in €{k) of type A. Then, by definition, 

^RC(A) = Catfc. (9) 

Ahfc 

These refined Catalan numbers have previously been studied by Haiman [ pj[ ] and Stanley pO[ ] in connec- 
tion with parking functions, and are known explicitly: 

^^(-'^ = (|A|-£(A) + l)!n.>i-.(A)!' ^''^ 

where rrii (A) is the multiplicity of i in A. Finally, given a pair of partitions A, /i h k, introduce the set of 
sequences of partitions 

«R(A,^)=={(A(^\...,A^('^))|A« h^„ A(i) U • • • U A(^('^» = A}. (11) 
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Thus sequences in 5H(A, /i) are obtained by breaking parts of /i in such a way that, after sorting, one 
obtains A. If CH(A, jjl) is non-empty, then A is said to be a refinement of /i. 

Theorem 1 Let tt G S{n) be a permutation of reducecF\cycle type fi \- k, and let A be another partition 
ofk. The number of primitive factorizations of?: of type A is 

J2 rc(a(i))...rc(aw^))). 

(A(i),...,A(*(a')))gK(A,ai) 

As an example, consider the permutation 

vr = (123456)(78910) £ S{10). 

This permutation has reduced cycle type fi — (5,3), so the length of a minimal primitive factorization of 
TT is eight. By Theorem y the number of minimal primitive factorizations of tt of type A = (3, 2, 2, 1) is 

RC(3, 2) RC(2, 1) + RC(2, 2, 1) RC(3) = 5 ■ 3 + 10 ■ 1 = 25. 

Theproof of Theorem |l|isbijective, and we refer the reader to our full-length article [ pOQ for details. As 
a corollary of this result, we obtain an elegant formula which counts the total number of minimal primitive 
factorizations of an arbitrary permutation. Gewurz and Merola's result is recovered as the case /i = (n) 
of this corollary. 

Corollary 1 Let -k G S{n) be a permutation of non-reduced cycle type /i h n. The total number of 
primitive factorizations ofir into n — i?(/i) transpositions is 

3 Primitive Factorizations and Jucys-Murphy Elements 

We will now give an algebraic explanation of the fact that the primitive factorization problem is central. 
Let C[S'(n)] denote the group algebra of the symmetric group, and Z{n) its center. The center is a 
commutative algebra with canonical basis {C^ : /i h n} consisting of the conjugacy classes of S{n)] for 
this reason we call Z{n) the class algebra. 

Let A denote the algebra of symmetric functions over C We define a specialization A -^ Z{n) as 
follows. For fc > 1, put 

Jfe := >, transpositions in S{k) — N^ transpositions in S{k — 1) = (1, A:) + • • • + (fc — 1, fc). (12) 

Thus Ji, . . . , J„ g C[S'(n)], with Ji = 0. The elements so defined are called Jucys-Murphy elements. 
They were introduced independently by Jucys [ |l^ ] and Murphy [ p2| ] . These simple elements have many 
remarkable properties, some of which we will make use of here. Diverse applications of the Jucys-Murphy 
elements are found in the work of Okounkov p3, Ga] and Okounkov and Vershik | 



''' Recall that the reduced cycle type of it is the partition obtained by subtracting one from the length of each of its cycles. Thus the 
size of the reduced cycle type of tt is the length of a minimal factorization of tt into transpositions. 
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Although { Ji, . . . , J„} % Z(n) for n > 3, the JM elements do belong to the Gelfand-Zetlin subal- 
gebra of C[5'(n)]. This is the maximal commutative subalgebra of C[5'(n)] generated by the class alge- 
bras Z{1), . . . , Z{n), where Z{1), . . . , Z{n — 1) are embedded in C[S{n)] in the canonical way. This 
is clear, since Jk is by definition the difference of an element of Z{k) and an element of Z{k — 1). 
Consequently, the JM elements commute with one another, and we may define the alphabet S„ = 
{{ Ji, . . . , J„, 0, 0, . . . }} and evaluate symmetric functions on this alphabet. It is a remarkable result 
of Jucys that, for any / G A, 

/(S„) e Z{n). (13) 

Thus symmetric functions of JM elements are central, and we have the JM specialization A — ;> Z{n). 

The JM specialization gives a very clean proof of the centrality of the primitive factorization problem. 
Let /ifc e A denote the complete homogoneous symmetric function of degree fc, i.e. 

hk^ ^ Xi^...Xi^, (14) 

the sum of all degree k monomials in the variables a;i, a;2, • • • • Then, by Jucys' result, /ife(S„) G Z{n) 
and we may write 

hki^n) = ^ ak,fj^Cfj, (15) 

fx\-n 

for some coefficients a^ p. On the other hand, we have that 

/ifc(s„) = 22 Jtf-Jtk 

2<ti<---<tk<n 



(16) 



2<ti<---<tk<nsi<ti s^Ktk 

= yj #{primitive factorizations of tt into k transpositions}?!. 

7reS(n) 

Thus, for any permutation n G S{n) of cycle type /i, we have that 

Qfc,/^ = 7^{primitive factorizations of tt into k transpositions}. (17) 

In other words, the combinatorial problem of counting primitive factorizations is equivalent to the alge- 
braic problem of resolving /ife(S„) with respect to the canonical basis of the class algebra Z{n). More 
generally, for any partition A one may consider the resolution 

mA(S„)=XfeA^^P, (18) 

where m\ is the monomial symmetric function of type A. We then have, for any tt G C^, 

bxi^ — #{primitive factorizations of n of type A}. (19) 



6 Sho Matsumoto and Jonathan Novak 

4 The Case of a Single Cycle 

The above algebraic encoding of the primitive factorization problem allows us to enumerate primitive 
factorizations of a full cycle tt S C(„) . This is thanks to the remarkable properties of JM elements in 
irreducible representations of C[S'(n)] , which amount to the fact that, while it is difficult to compute the 
the coordinates of /(S„) with respect to the class basis of -Z(n), it is easy to compute its coordinates with 
respect to the character basis of Z{n). 

The following remarkable result is due to Jucys [|l8|], see [Q] for a complete proof. Given a Young 
diagram A and a cell D in A, recall that the content of A is defined to be the column index of A less its row 
index. Let us associate to each Young diagram A its alphabet 

A^ = {{c(n) : D e A}} (20) 

of contents. Let H\ denote the product of all hook-lengths of A. Let {x"^ : A h n} be the characters of 
the irreducible representations of S{n), which form a basis of Z{n). Then, for any symmetric function 
/ G A, the character expansion of f[^n) is 



X\-n 



/(s„) = V :^^x^ (21) 



Thus, at the combinatorial level, the character expansion of f{^„) is implemented by the substitution rule 
S„ -^ Ax. 

Let us use the above character expansion result to enumerate primitive factorizations of a full cycle 
TT e C(„) of any given length. We akeady know that the number of minimal primitive factorizations of 
(12 . . .n), i.e. those consisting of n — 1 transpositions, is the Catalan number Cat„_i . We will now 
solve the problem when the number of transpositions used is allowed to be arbitrary. As we will see 
momentarily, this problem has an unexpectedly simple and beautiful solution. 

Let z be an indeterminate, and form the generating function 

$(z;n) = ^/lfe(S„)z^ (22) 

fe>0 

This generating function is an element of the algebra Z(n)[[z]] of single-variable formal power series with 
coefficients in the class algebra Z{n). By Jucys' character expansion result, we have 



*(-;-) = E E 



hk{Ax) A j^ 



fc>0 "- \\-n ^ 



-X z 



= T.(T.''^i^^>')jrx (23) 






E- 



^_^HxUnexi^-<^)^) 
where the last line follows from the generating function 



E 

fe>0 



/ife(xi,...,x„)z'= = TT- (24) 

r"--"- 1 — XiZ 
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for the elementary symmetric functions. Note that this computation shows that the generating function 
$(z; n) is actually a rational function over Z{n). 

Now, given /i h n, in order to obtain the generating function 

^MW-II«fe,M^'' (25) 

we simply take the corresponding traces of the conjugacy class C^ in each irreducible representation, to 
obtain a rational function over C : 

Up until this point, the partition /i h n has been generic, but now we restrict to the special case /i — (n), 
the partition of n with a single part. A classical result from representation theory informs us that the trace 
of C(„) in an irreducible representation can only be non-zero in "hook" representations: 

^A(^^)^ f(-ir, ifA = (n-r,r) _ ^^7) 

I 0, otherwise 

Now, the content alphabet of a hook diagram may be obtained immediately, 

Ain-r,in = {0, 1, . . . , n - r - 1} U {-1, . . . , -r}. (28) 

so that 

^(n){z)^y^ T • (29) 

For example, if n = 4, this is a rational function of the form 



const. const. 

^^'>^^' "(l-z)(l-2z)(l-3z) ^ (l-z)(l-2z)(l + z) 
const. const. 



(30) 



(1 -z)il + z)(l + 2z) (1 + z)il + 2z)(l + 3z) ■ 
Thus, as an irreducible rational function, $(„) (z) has the form 

'^(")W = T?^^7r^ ^31) 

n,=i (1 - «^2^) 

where cq, . . . , c„_i G C are some constants to be determined momentarily. 
Before finding the above coefficients, let us consider the generating function 

^ Y,hg{l\...,n')ua. (32) 
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The coefficients in this generating function are complete symmetric functions evaluated on the alphabet 
{1^, . . . ,n^} of square integers. Reason dictates that they ought to be close relatives of the Stirling 
numbers 

S{n + g,n) = hg{l,...,n). (33) 

The Stirling number S{a, b) has the following combinatorial interpretation: it counts the number of parti- 
tions 

{l,...,a} = ViU---UVb (34) 

of an o-element set into b disjoint non-empty subsets. Stirling numbers are given by the explicit formula 



T{n + g,n) = hg{l^,...,n^) (36) 



The numbers 

are known as central facto rial numbers. The central factorial numbers were studied classically by Carlitz 
and Riordan, see [BU Exercise 5.8] for references. They have the following combinatorial interpretation: 
T{a, b) counts the number of partitions 

{l,l',...,a,a'} = T/iU---UVb (37) 

of a set of a marked and a unmarked points into b disjoint non-empty subsets such thatpj, for each block 
Vj,\fi is the least integer such that either i or i' appears in Vj , then {i,i'} C Vj. Central factorial numbers 
are given by the explicit formula 

T'-M^^p-.r'j.^^^, (38) 

Now let us determine the unknown constants cq, . . . ,c„_i. By the above discussion, the generating 
function $(„) (z) has the form 

$(„) (z) = (CO + ciz + • . • + c„_iz"-i) J2 T{n -l+g,n- l)^^^. (39) 

3>0 

On the other hand, by the results of the previous sections, 

k>0 

= 2_] On- i+2g, (n) z"^ ^+^3 , siucc cvcry permutation is either even or odd, (40) 

= Cat„_i z""^ + a„+i^(„)z"+^ + . . . , since a„_i^(„) = Cat„_i . 

Consequently, we must have cq = ■ • ■ = Cn-2 = 0, c„_i = Cat„_i, and we have proved the following 
result. 

'"' Balint Virag gave a colourful description of this condition, which is actually quite a useful mnemonic; "the most important guy 
gets to bring his wife." 
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Theorem 2 For any g > 0, the number of primitive factorizations of (12 . . . n) G S{n) into n — 1 + 2g 
transpositions is 

Cat„_i ■T{n — 1 + 5,ri — 1), 

where T{a, b) denotes the Carlitz-Riordan central factorial number Equivalently, we have the generating 
function 

/ N _ ^at„— 1 z 

^"'^^^ " (l-l2z2)...(l-(n- 1)2^2) ■ 

5 Primitive Factorizations and Matrix Models 

Finally, we come to what is perhaps the most intriguing aspect of the primitive factorization problem: its 
connection with matrix models. The theory of matrix models has its origins in an area of mathematical 
physics known as quantum field theory, see [g|] for a solid introduction. For our purposes, the following 
grossly oversimplified description of matrix model theory suffices: 

1 . Pick an interesting subset S{N) of the space M (N) of all A^ x A^ complex matrices. 

2. Put an interesting probability measure r/ on S{N). 

3. Select an interesting random variable (measurable function) / : S{N) — > C. 

4. Compute the expected value (/) of the random variable / (possibly after rescaling) as a power 



series m -^ 



eM) 



if)^ I J[MUdM)^^^. (41) 



5. Realize that the coefficients Eg (/) occurring in the above perturbative expansion have an interesting 
combinatorial interpretation. 

This informal discussion is meant to convey the impression that, from a combinatorial perspective, matrix 
integrals may sometimes play the role of generating functions. It sometimes happens that a traditional 
generating function is difficult to obtain, but that by running the above steps in reverse one can concoct 
a matrix integral which encodes a sequence of interest. Furthermore, some matrix models have special 
features that are very useful, and these features can be used to extract combinatorial information in the 
generating function spirit. For example, it might be that an integral of interest can be exactly evaluated, 
thereby yielding an explicit generating function for the sequence it encodes (see [ p^ for a famous exam- 
ple). Even if this is not the case, it often happens that matrix integrals interact well with more advanced 
analytical tools, such as orthogonal polynomials or integrable systems of differential equations (see e.g. 
[JSJ]), and this may again yield insight into combinatorial structure. 

Let us present a matrix model for primitive factorizations. For our space of matrices we select the 
group U{N) of N X N complex unitary matrices. Since U{N) is compact, it carries a unique left and 
right translation invariant probability measure, the Haar measure dU, which we take for our probability 
measure of interest. Now we will select an interesting random variable, or rather class of random variables, 
/ : U{N) -^ C. It would certainly be nice if we could compute the expected value {P{U, U~^)) of 
any polynomial function of the entries of U and U^^ — U*, since we can approximate a large class 
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of functions on U{N) by polynomials in matrix coefficients (think Stone- Weierstrass/Peter-Weyl). By 
linearity of the integral, it suffices to consider the case where P is a monomial. Furthermore, an easy 
argument using the invariance of the Haar measure shows that the expected value of such a monomial will 
be zero unless P is of equal degree in the entries of U and U* (think integrals of the form / z'"z"dz over 
the unit circle U{1)). Thus we need only consider integrals of the form 

(Wi(l)i(l)"i'(l)j'(l) • • • Wi(n)j(n)Ui'(n)i'(n)) = / "i(l)j(l)"i'(l)j'(l) • • ■ Ut{n)j(n)Ui' {n)j' (n)dU, (42) 

JU(N) 

where the lowercase u^'s are matrix elements and i, j, i' , j' : {1, . . . ,n} —!■ {1, . . . , N} are functions. 
Integrals of this form are called n-point correlation functions of matrix elements. They are actually of 
considerable interest in mathematical physics [|| |l2[ |T]] and free probability theory [g|. It is known 
that, provided N > n, the computation of the n-point functions can be reduced to the computation of 
"permutation correlators" 

(uilUi^d) ■ ■ ■ UnnUn^{n)), TT G S{n). (43) 

Finally, recall that we have defined 

afc,Ai — #{primitive factorizations of vr into k transpositions}, (44) 

and that this quantity can be non-zero only for k of the form k — n — £{fj,) + 2g. Therefore let us introduce 
the notation 

<ig,p, '■= 0-n-i{ii)+2g,ii- (45) 

It is not unreasonable to think of Og ^ as a combinatorially motivated analogue of the usual Hurwitz 
number [B] hg^^, obtained by replacing the transitivity constraint with the primitivity constraint. 

Theorem 3 Let ^ be a partition of n and let vr £ C^ be a permutation of cycle type fi. Then, for any 

N >n, 

S>0 

Unfortunately, we will not be able to say much about the proof of this result here. Suffice to say that 
that Theorem ^ arises from two points of view regarding the orthogonal projection of M (A^)®" onto the 
commutant 

CuiN){n) ^{Te 7W(iV)®" : U®''T = TU®"" VC/ G U{N)]. (46) 

The first point of view involves the permutation correlators ( p3| ) and the second involves understanding 
the element [N + Ji)^^ . . . {N + J„)~^ in the left-regular representation of C[5(n)]; the equivalence 
of the two is, in a sense, a manifestation of the Schur-Weyl duality between the representation theories 
of S{n) and U{N). We refer the interested reader to our articles [20, ESJ for further details regarding the 
proof of Theorem bL and its applications. 
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Let us conclude with the following comparison of unrestricted and primitive factorizations of a full 
cycle: 






(2.g)! 



(2g)! 






The first of these formulas is due to Jackson [ p^ (see also [p8|]), while the second is a consequence of 
Theorem together with Riordan's exponential generating function for the central factorial numbers. 
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